January,2001 



O 
O 

(N 

C 

c3 



> 



6 

c3 



13 



Toy model for a two-dimensional accretion disk 
dominated by Poynting flux 

Hyun Kyu Lee* 

Department of Physics, Hanyang University, Seoul, 133-791, Korea 

and 

Asia Pacific Center for Theoretical Physics, Seoul 130-012, Korea 

Abstract 

We discuss the effect of the Poynting flow on the magnetically dominated thin 

accretion disk, which is simplified to a two-dimensional disk on the equatorial 

plane. It is shown in the relativistic formulation that the Poynting flux by 

the rotating magnetic field with Keplerian angular velocity can balance the 

energy and angular momentum conservation of a steady accretion flow. 
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I. INTRODUCTION 



The Poynting flux model |T| |2|] suggested for an accretion disk with the ordered magnetic 
fields has been considered to be one of the viable models for the astrophysical jets || 0|. In 
contrast to the hydrodynamic jets in which the energy and angular momentum are carried 
by the kinetic flux of the matter, the Poynting flux is characterized by the outflow of energy 
and angular momentum carried predominantly by the electromagnetic field. 

In the non-relativistic formulation, Blandford PJ suggested an axisymmetric and steady 
solution for the Poynting outflow assuming a force-free magnetosphere surrounding an ac- 
cretion disk. The poloidal field configuration for a black hole in a force- free magnetosphere 
has been discussed recently by Ghosh || in the relaticistic formulation, in which the pos- 
sible forms of the poloidal magnetic fields are suggested. The developments of the ordered 
magnetic field and the Poynting outflow on the disk have been studied by many authors 
H |7j H and recently Ustyugova et al. jp perform a axisymmetric magnetohydrodynamical 
simulation to show that the quasi-stationary and approximately force-free Poynting jet from 
the inner part of the accretion disk is possible. 

The Poynting flux in a system of black hole-accretion disk recently has also been studied 
in connection to the gamma ray bursts || [Kj |TI |. One of the main reasons is that the 
Poynting flux carries very small baryonic component, which is essential for powering the 



gamma ray bursts ||12|| . The evolution of the system is found to be largely depend on the 
Poynting flux on the disk [^4[ . The relativistic effects on the accreting flows close to a 
black hole has been discussed by Lasota [|T5 1 , Abramowicz et al. fllfj and recently by Gammie 



and Popham [|L7] in which the relativistic effects on the slim-disk with the vertical structures 
averaged are discussed in detail with viscous stresses. However the effect of relativity on the 
accreting flow dominated by the Poynting flux has not been discussed well in depth so far. 

The purpose of this work is to study the effect of the Poynting flux on the accretion 
flow in the relativistic formulation. In this work, we consider a toy model for a magnetically 
dominated thin accretion disk, which is assumed to be a two-dimensional disk located on the 



2 



equatorial plane with a black hole at the center. To see the magnetic effect transparently it 
is also assumed that there is no viscous stress tensor in the disk and there is no radiative 
transfer from the disk. We develop a relativistic description of the two-dimensional model 
for a magnetically dominated thin accretion disk in the background Kerr geometry. It is 
shown that the energy balance of the accretion disk for a stationary accretion flow can be 
maintained by the Poynting flux, provided that the poloidal magnetic field is rotating with 
the same Keplerian angular velocity Qk on the disk. 



The stress-energy tensor Tfif of the matter in the disk can be given in the general form 



where p m , IT, and p are the rest-mass density, internal energy and pressure respectively. S^ v 



In this work we simplify the accretion disk to be non-viscous^^ = 0), cool(p = 0, II = 0) 
and non-radiative (g M = 0) disk to investigate the effect of the magnetic field transparently. 
It is also assumed that there is negligible mass flow in the direction perpendicular to the 
disk: u — 0. Then the stress-energy tensor is given simply by 



II. TWO-DIMENSIONAL ACCRETION FLOW 



T% = (p m + p + Il)uV + pg^ + + u»q v + u v q\ 



(1) 



is the viscous tensor and g M is the radiative energy flux [[UJ . 




p m u^u u . 



(2) 



The energy flux is obtained by 



(m) 



-p m u u^, 



(3) 



and the angular momentum flux by 



(m) 



p m u^ur. 



(4) 



For an idealized thin disk on two-dimensional plane we take 
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Pm = ^{9 - 7T/2), (5) 

where a m is the surface rest-mass density. In this work, the background geometry is assumed 
to be Kerr metric(see Appendix A) with rotating black hole at the center and we adopt the 
natural unit G = c = 1 . 

The rate of the rest-mass flow crossing the circle of radius r is given by 

/ ap m u r ^-j=d6d(f) = 2ir<j m pu r , (6) 

Jr V A 

which defines the mass accretion rate M + by 

M+ = -27ia m pu r . (7) 



This expression is identical to that derived in [|16[] and [17], in which the vertical structures 
are averaged. The stationary accretion flow implies that M + is r-independent. 
Using eq.(|) and (f|) the radial flow of the matter energy at r can be given by 

2 

/ aE\ m y-^=dQd(\) = -lTxo m pu r u = u M+ (8) 

Jr V A 

and the radial flow of the matter angular momentum by 

2 

/ a£ r , m) ?-^Ld6d<j) = 1i\o m pu r = -u^ M + . (9) 

Jr \ A 



X A 

For the magnetosphere outside the accretion disk, there are electromagnetic currents, 
J M , which can flow through the central object at the center and also along the magnetic 
field lines anchored on the disk. In this work we suppose a situation in which the currents 
flow into the inner edge of the disk through the central object and the currents flows out 



along the magnetic field lines from the disk as discussed by Blandford and Znajek [g] |]T 
The continuity of the currents(current conservation) necessarily requires surface currents in 
radial direction on the disk. Similar consideration has been applied to the black hole horizon 



as discussed by Thorne et al. [ |H|j . 

In general the structure of the electromagnetic field with discontinuity at 9 = tt /2 plane 
can be reproduced by assigning the surface charge and current on the plane in addition to 



the bulk charge and current distributions which terminate on the disk. The surface charge 
density a e and surface current density K % can be defined systematically using the procedure 
suggested by Damour [^] (See Appendix B for details): 

a e = --, (10) 
K f = -—B\ K* = —B f . (11) 

III. POYNTING FLUX 

The energy and angular momentum of the disk can be carried out by Poynting flux along 
the magnetic field lines anchored on the disk. In our simple model, Poynting flux is the main 
driving force for the accretion flow. We calculate the Poynting flux on a two-dimensional 
disk assuming a force-free magnetosphere around the accretion disk. 

Using the Killing vector in t-direction, 

e = (1, 0, 0, 0), (12) 
we can define the energy flux £ M from the energy momentum tensor, T pv 

£ » = = (a 2 - w 2 (3)T p0 - w 2 (3T^, (13) 

= ( 14 ) 

where 

T" u = ^-(F p p F up - -g^F^F" ). (15) 

Since we are interested in ^-direction (normal to the disk on the equatorial plane ), we 
consider S e for which the second term in eq.(|T^) vanishes: 

S e = —{a 2 - m 2 (3)F o F op - —w 2 (3F e F^. (16) 
An p 4tx p 

Using the identities given in Appendix A, it can be rewritten in terms of electric and magnetic 
fields as given by 



' [aE x B\ e + pw(E § E* + B^B § )\. (17) 



The power measured at infinity can be obtained using eq.(Q) by 



PLr 9y = J o£° dZ e = J ^[a(E x B) B + w(3(E § E* + B^B^^Ldrd^. (18) 



--n/2 J 4vr VA 

When the tangential components of the electromagnetic field are multiplied by the laps 



function a as in the membrane paradigm flTJ], we can hide a in the integrand. 



For the steady and axisymmetric case which we are interested in , E^ = 0, we get 

£ e = —(-aE f B^ + f3wB^B § ). (19) 



The first term in the integrand in eq.(|I~8"D can be rewritten in terms of the surface current 
density: Using eq.flTTD we get 

— {E x B) e = E ■ K. (20) 
An 

For the current in the direction of tangential electric field on the disk there is a energy 



dissipation into the disk surface. This is what one can expect on the black hole horizon |L9 
However for the current in the opposite direction this term corresponds to the electro-motive 
force and it is the case for the accretion disk on the equatorial plane discussed in Blandford 
and Znajek |2| |l£|as well as in this work. 

Similarly the third term in eq.(^)(equivalently the second term in eq. (|T9|)) can be rewrit- 
ten in terms of the surface current using eq.([7T|) as given by 

—w(3B*B § = wuK f B § . (21) 

We can see that this is the magnetic braking power on the rotating body with an angular 
velocity u = —f3. 

Using the Killing vector in ^-direction for the axial symmetric case, 

= (0, 0, 0, 1), (22) 
we can define the angular momentum flux £ M from the energy momentum tensor, T^ v , 



£M = T ^ Vu = W ^T"° + w 2 T^, (23) 
;„ = 0, (24) 

Then we get the flux in ^-direction (normal to the disk on the equatorial plane ), C , 

C? = — —B®B^ = -wK f B\ (25) 
4?rp p v ' 



which is nothing but a torque applied on the surface current density K r . Eq . (|17D can be 
written as 

£ e = — aE x B\ g -f3C e . (26) 

4:TCp 

In this work we assume that there is a sufficient ambient plasma around the disk to 
maintain the force-free magnetosphere @ Q . The force-free condition for the magnetosphere 
with the current density J M is given by 

V = 0. (27) 
The magnetic flux, through a circuit encircling <fi = — > 2ir, 

* = £ A^dcfi = 2nA h (28) 

defines a magnetic surface on which A^r, 9) is constant and therefore is characterized by the 
magnetic flux \l/ contained inside it. From the force- free condition it can be shown that Aq 
is also constant along the magnetic field lines and the electric field is always perpendicular 
to the magnetic surface. We can define a function flp( r , 9), 

dA = -SlpdAfr (29) 

which is also constant along the magnetic surface. VLp can be identified as an angular velocity 



of magnetic field line on a magnetic surface |L8] [[19 
Then we get from the force-free condition 



E = — — {VLf + 0)<j> x B, (30) 
a 



and the Poynting flux perpendicular to the disk can be written as 

£ " = -i^ B " B * = a "^y (31) 



as in fl8fl . 



IV. ENERGY AND ANGULAR MOMENTUM CONSERVATION 

In the idealized thin disk considered in this work, it is assumed that there is no radiative 
transfer or viscous interaction which balance the radial flow of the energy and angular 
momentum of the accreting matter. Hence the electromagnetic field anchored on the disk 
is responsible for the conservation of the total energy and angular momentum. Since we 
assume the idealized two-dimensional disk, there is no radial flow of electromagnetic energy 
and angular momentum unless there is any singular structure of electromagnetic field on the 
disk. 

Let us consider a circular strip at r with infinitesimally small width Sr. Using eq.(p5[) 
the angular momentum flux of electromagnetic field into the ^-direction is given by 

rr+Sr 

AL = 2 aC B <nie = -B e B^pw5r, (32) 



where the factor 2 is introduced to take account of both sides of the disk. It is balanced by 
the change of the angular momentum of matter given by 

AL m = -u^M + \ r r +Sr -> -^M + 5r. (33) 

Then we get 

- B § B^>pw 
+ du^/dr 

In the non-relativistic limit can be identified as non-relativistic Keplerian angular mo- 



mentum: Vt K = J Mj r 3 and we get 



B e B*r 2 B e B* , x 
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which is identical to that used in [H31 and 



Similarly the energy flux of electromagnetic field in ^-direction is given by 

AE = 2 f +5r aS e dY> e = —VtpB®B^pw5r (36) 



= tt F AL, (37) 

where factor 2 is introduced for the same reason as in the angular momentum flux. It is 
balanced by the change of the energy of the matter given by 

AE m = -(- Uo )M + \ r r +dr -> ^M + 6r, (38) 

which is possible with Qp given by 

n, = (39) 

du^/dr 

where eq.(p4|) has been used. 

It is not a trivial task to find the solution of the full relativistic MHD equations satisfying 
eq.(|7|) and ( p4|) which leads to 

B 6 B* = -2na m pu rd ^, (40) 
dr 

together with eq. (]3"5|) . Without solving the full MHD equation, the magnetic fields and its 
angular velocity Vtp can be considered to be unknown parameters as well as ttj. 

As a first trial, we assume uq and to be those of the stable orbit of a test particle 



around a Kerr black hole 21 



r 2 - 2Mr + ay/Mr 
u = — , 



r\ir 2 - 3Mr + 2a\[Wr 



y/Mr(r 2 - 2aVMr + a 2 ) , . 

«0 = , - , (41) 

r Vr 2 - 3Mr + 2aVMr 



although it is natural to expect non-negligible effect of the magnetic field on the stable 
orbit. The straightforward calculation, eq.(^), shows a very interesting result that the 
angular velocity Qp of the magnetic field is the same as the Keplerian angular velocity 



n F = n K , (42) 



where 



n K =J- t== m (43) 

V r 3 i + a^/M/r 3 

In fact one can derive eq.(|4*2]) formally without using the explicit form of eq. (STj) , since it 
can be shown that 

u0 duo + A = Q (44) 
ar dr 

holds in general for a stable circular orbit on the equatorial plane. This shows that the 
energy balance can be realized by the magnetic field rotating rigidly with Keplerian angular 
velocity in this simplified two-dimensional disk model. 

It is also possible to show that eq.(4.2) in ||| can be obtained from eq. (|34]) in the non- 
relativistic limit. Hence for a Newtonian limit or for the outer radius of the thin disk, the 
configuration of the electromagneic field suggested by Blandford [Q is consistent with our 
result when the the magnetic field angular velocity Q F is given by VLk- 

V. POWER OUT OF A MAGNETICALLY DOMINATED THIN ACCRETION 

DISK 

The total power out of the disk, Pdisk, can be calculated by integrating eq. (|36|) from the 
inner most stable point(r iri ) to outer most edge of the disk with Poynting flux(r OU i), 



disk 



-n F B e B' t, pw)dr, (45) 



where — sign is referring the outward direction on the disk. 

Using the accretion rate eq. (|34T) , which is independent of r for a stationary accretion 
flow, and the energy balance condition eq.(^), we get 

^ f r ° ut ^ ,v dvbfk , • r r ™t d(—u ) , , . 

P disk = / n F M + —±dr = M + / -^^dr, (46) 
Jr in ur Jr in ctr 

= [-u (r out )]M + - [-u {r in ))M + . (47) 
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The first term corresponds to the energy rate into the disk at the outer most edge of the 
disk and the second term corresponds to the energy accretion rate at the inner edge of the 
disk to the central object. This is the expected result from the energy conservation. 

Eq.pTD indicates that the power out of the disk in this toy model is not dependent 
strongly on the details of the electromagnetic field configuration. The power can be cal- 
culated once M + is known at one point r on the disk. For example, suppose that the 
exact solution of the relativistic equation can be approximated at a large distance r$ by the 
configuration suggested by Blandford @, then we can use the relation, for example, 

= 2Q F rB § , (48) 

to get the accretion rate expressed in terms of the magnetic field component perpendicular 
to the disk, 

M + = A(B § (r )r ) 2 . (49) 
Then the total power for a disk with Poynting flux extended to r out = oo can be given by 

P disk = 4(1 - [-u (r in )})(B 8 (r )r o y, (50) 
where we take — -u (oo) = 1. 

VI. DISCUSSION 

In this work, we discuss a toy model for the magnetically dominated thin accretion disk. 
Assuming a two-dimensional disk dominated by the Poynting flux where the viscous stress 
and radiative transfer are ignored, the accretion flow in two-dimension is discussed in the 
background of Kerr geometry. We have demonstrated that the stationary accretion can be 
possible with the co-rotating magnetic field with the same Keplerian angular velocity as that 
of matters in the disk in the stable orbit. Also it is observed that the solution proposed by 
Blandford [[J is consistent with our result in the non- relativistic limit ( or at a large distance 
r from the center with sufficiently small a/r and a/M). It is shown in this toy model that 
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the total Poynting power out of the disk depends on the accretion rate and the energy at 
the inner most stable orbit but not on the details of the electromagnetic field configuration. 

Related issues to be discussed in the future are the possible solutions of magnetic field 
configuration(not only the poloidal component of the magnetic field discussed by Ghosh || 
but also the toroidal component) analogous to in the non-relativistic limit and the effect 
of the magnetic field on the energy(— uq) and angular momentum)^) of a particle in the 
stable orbit which has not been taken into account in this work. 

We have thus far discussed only one specific aspect, Poynting flux, of the accretion 
disk. For this two-dimensional disk model to be realistic and viable one, there are many 
obstacles to overcome. For example the viscous stress and radiation transfer should be 
included, which naturally leads to the questions on the steady state || and the thin disk 
approximation assumed in this work. Moreover recent works |22|] [^3| J24| on the accretion 



flow toward a black hole as a central object and on the magnetic coupling |14[ between the 
disk and the black hole seem to indicate that the physics is more complex than the simplified 
two-dimensional model discussed in this work particularly near the inner region of the disk. 
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Appendix A : Kerr geometry 



In this Appendix, the useful identities expressed in terms of explicit orthonormal com- 
ponents defined by ZAMO are listed. Using the Boyer-Lindquist coordinates |S5| in the 
natural unit G = c = 1, the metric tensor for Kerr Geometry f26|, g^, is given by 
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(-(a 2 -w 2 (3 2 ) w 2 (3\ 






P 2 

A 














P 2 





zu 2 [3 








w 2 



where 



a = — — — , p = — -, 00 = —smt), 

£ 9H> P 



A = r 2 + a 2 -2Mr, p 2 = r 2 + a 2 cos 2 6, £ 2 = (r 2 + a 2 ) 2 - a 2 A sin 2 fl, 



(51) 
(52) 



ZAMO's four velocity(a time-like unit vector orthogonal to the t-constant surface: dx a U a = 
0) is given by 



U„ = -a(l, 0, 0, 0). 



(53) 



The electromagnetic field tensor can be expressed by the electric and magnetic fields as 
given by 



\=B*, F e = \(apE § + PpwB f ) 
A p 2 



F 0r = ^ 

ap 



aw ' 



ap 



zup 



(54) 
(55) 



Appendix B : Surface current on a two-dimensional accretion Disk 

Consider a conserved current J 11 defined by 

jv. = jn Y (6 - tt/2) + f, (56) 
where the conserved current J M is the bulk current density satisfying the Maxwell equation, 

F^ = 4ttJ», J^ = (57) 

and Y(9 — ir/2) satisfies 

Y{6 - tt/2),,, = -5 e ,5(9 - tt/2) (58) 
13 



From the conservation of the current, 



J% = 0, (59) 

we get, 

^ = ^;^-^/2)- (60) 

Then we can identify 

f = 1.^6(6-^/2), (61) 



since F vanishes identically. It is an analogous expression as in |20 |. 
Now the charge density defined by 

p e = aj° = a— 5(6-n/2), (62) 

can be rewritten as surface charge density a e : 

J p e dV = J a e —j=zudrd(f). (63) 



A 

Then we get the surface charge density in terms of the electric field given in Appendix A: 

°e = ^ (64) 

It is Gauss' law on the disk plane: 

E § = -Ana e . (65) 
Similarly using the current density defined by 

~f=f-f(0, 0,-/3), i = r, 6, <\> (66) 

we get 

> = y —6{6-rr/2), ~f = -—5(6 - n/2), (67) 



4/T p 2 4tt pw 
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where j e = by construction. The radial surface current density K r and the surface current 
density in 0— direction on the disk can be defined by 

j afdUr = ~j- j aB *™ d( t> = J ®K f wd(f), (68) 

to get 

JC = --B*, = —B f , (70) 

An 4ir 

where dT^i is the corresponding surface element. 

This result can be summarized by Ampere's law on the disk surface: 

B = -AnK x 0. (71) 



15 



REFERENCES 

[1] R.V.E. Lovelace, Nature 262, 649(1976) 
[2] R.D. Blandford, MNRAS 176, 465(1976) 

[3] M.C. Begelman, R.D. Blandford and M.J. Rees, Rev. Mod. Phys. 56, 225(1984) 

[4] G.V. Ustyugova, R.V.E. Lovelace, M.M. Romanova, H. Li and S.A. Colgate, ApJ 541, 
L21(2000) 

[5] P. Ghosh, MNRAS 315, 89(2000) 

[6] P. Ghosh and M.A. Abramowicz, MNRAS 292, 887(1997) 
[7] M. Livio, G.I. Ogilvie and J.E. Pringle, ApJ 512, 100(1999) 
[8] C. Fendt, Astron. Astrophys. 319, 1025(1997) 

[9] H.K. Lee, R.A.M.J. Wijers and G.E. Brown, Phys. Rep. 325, 83(2000) 
[10] H.K. Lee, G.E. Brown and R.A.M.J. Wijers, ApJ 536, 416(2000) 
[11] L.-X. Li, Phys. Rev. D61, 084016(2000) 
[12] T. Piran, Phys. Rep. 314, 575(1998) 

[13] H.K. Lee and H.K. Kim, J. Korean Phys. Soc. 36, 188(2000) 
[14] L.-X. Li, B. Paczynski, ApJ 534, L197(2000) 

[15] J. -P. Lasota, in "Theory of Accretion Disks 2" ed. W.J. Duschl et al.(Dordrecht:Kluwer, 
1994) 

[16] M.A. Abramowicz, X.-M. Chen, M. Granath and J.-P. Lasota, ApJ 471, 762(1996) 

[17] C.F. Gammie and R. Popham, ApJ 498, 313(1998) 

[18] R.D. Blandford and Z.L. Znajek, MNRAS 179, 433(1977) 

16 



[19] K.S. Thorne, R.H. Price, and D.A. Macdonald, "Black Holes; The Membrane Paradigm" 
(Yale University Press, New Haven and London, 1986) 

[20] T. Damour, Phys. Rev. D 18, 3598(1978) 

[21] S. Shapiro and S.A. Teukolsky, "Black Holes, White Dwarfs, and Neutron Stars" (John 
Wiley & Sons, New York, 1983) 

[22] C.F. Gammie, ApJ 522, L57(1999) 

[23] J.H. Krolik, ApJ 515, L73(1999) 

[24] B. Punsly, 1998 ApJ 506, 790(1998) 

[25] R.H. Boyer and R.W. Lindquist, J. Math. Phys. 8, 265(1967) 
[26] R.P Kerr, Phys. Rev. Lett. 11, 237(1963) 



17 



